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The standard model is a successful but limited theory. There is significant theoretical moti- 
■ vation to believe that new physics may appear at the energy scale of a few TeV, the lower end 



of which is currently probed by the Fermilab Tevatron Collider. The methods used to search 
for physics beyond the standard model in a model independent way and the results of these 
searches based on 1.0 fb~^ of data collected with the D0 detector and 2.0 fb~^ at the CDF 



^ ' detector are presented. 

^T) • 1 Introduction and Strategies 

' Particle physics is at a stage where there is no unique way forward. The standard model of 

■ particle physics has been remarkably successful: all fundamental particles predicted by this 

model have been discovered, with the notable exception of the Higgs boson. Despite its success, 
there are strong motivations from the theory to expect new physics at energies at or just above 
^ ' the electroweak scale. 

Assuming that new physics does exist, we do not know what this new physics is and thus 
precisely how to search for it. There are many theories and models which predict differences 
between reality and the standard model. But generally these theories and models do not give 
precise energy and phase space regions to search for new physics. Motivated by this, a scan 
over many channels was performed at each experiment to look for significant deviations from 
the standard model. The CDF experiment searches over all high-pT data. At DO, it was found 
that the background model was most developed for the case of final states containing lep tqns. 
Similar approaches to search for new physics have been applied to data from D0 Run I^tHSEl 
and HI at heraSI. 

At CDF and D0 , several methods are used to search across large portions of the data set. 
First, exclusive final states are tested for agreement in event counts using the vista algorithm. 
VISTA also searches many histograms for each of the final states for disagreement. In sleuth, 
one variable, {T,pT,obj) + MET, is used for the search, focusing on the tails of the distributions. 
Finally, at CDF, masses of all object combinations are searched for discrepancies in bump 

HUNTER. 

focuses on consistency and reach. The objects used in all final states have the same 
set of object cuts. The standard model background implementation comes from Monte Carlo, 



primarily PYTHi^ and MadEveN':^^. After the Monte Carlo is pushed through the detector 
simulation, additional corrections to fix modeling issues are performed in a standard, rigorous 
way. All of the final states and histograms are put into one immense fit, using as few parameters 
as possible to provide basic agreement across final states. The fit uses all of the exclusive final 
states defined in vista with the histograms provided for those final states. In practice, it is 
often one histogram that dominates the value found for any of the individual parameters. CDF 
introduced 43 correction parameters to minimize final state discrepancies. The parameters are 
further constrained by adding information external to that used in the vista high-pj^ data 
sample. The correction factors obtained from the fit are found to be consistent with those 
derived in control regions of other analyses at CDF. 

As at CDF, D0 runs Monte Carlo simulations to predict the standard model expectation. 
D0 uses mostly ALPGEiPand pythia Monte Carlo. D0 models multijet background from data 
rather than Monte Carlo. To estimate these backgrounds, some of the object selection cuts 
are reversed in data to produce samples with electron, muon, or tau objects which are mostly 
misidentified jets from multijet backgrounds. 

D0 then divides the whole data set into seven nonoverlapping final states. In each of these, 
previously determined standard weights are applied using well-understood areas of phase space 
dominated by particular standard model processes to account for necessary corrections to the 
Monte Carlo and detector modeling. These seven states are defined by their lepton content 
and are inclusive in jets and additional objects. A fit is then performed for each of these 
states to obtain the scale factors which reproduce the distributions of the selected data with the 
background from Monte Carlo and multijet background determined from data. 



2 General Searches Using vista 

VISTA performs two checks, a normalization-only check on the number of events in each exclusive 
state, and a shape-only analysis of histograms within a state by calculating a Kolmogorov- 
Smirnov statistic (and resulting fit probability). Both of these numbers require additional 
interpretation, because of the number of trials involved. When observing many final states, 
some disagreements are expected merely due to statistical fiuctuations in the data. Therefore, 
probabilities calculated with both of these methods are corrected to reflect this multiple testing. 

At CDF, there are no discrepancies in event counts over the 399 states checked. Additionally, 
555 1-D shape histograms show discrepancy after trials. Many of these can be attributed to a 
Monte Carlo "3-jet effect". The "3-jet effect" is seen in several types of histograms across 
many of the multijet final states. These include the ATZ{j2, js) and the jet mass distributions. 
Examining different parton shower generators show that predictions for this region of phase 
space are not consistent. These predictions have not yet been compared against LEP I data 
to provide reasonable constraints on this effect. The distributions of final state counts and 



histograms are shown is Figs. 1(a) and 1(b) 



At D0 , of the 180 distributions, four show significant event count discrepancies. These are 
the final states + 2 jets + I^t with a converted probability of 9.3a after trials correction, /x 
+ 7 -|- 1 jet -|- I^T with 6.6(T, -|- I^t with a discrepancy of 4.4(t and ^'^ ^~ + 7 at 4.1(7. 

Additionally, 24 histograms show shape discrepancies. The distributions of final state counts 



and shape discrepancies can be seen in Figs. 1(c) and 1(d) 



Two of these states are directly related to oversimplified modeling of the photon misidentifi- 
cation rate. The ^ + 2 jets + IpT final state discrepancy shows an excess of events with a muon 
at r/ > 1.0. The excess points to an oversimplification in our approach to trigger efficiencies. 
The proportion of events that are brought in by single muon vs. muon plus jets triggers changes 
significantly as we increase jet multiplicity. These triggers introduce r/-dependent efficiencies 
which are not fully incorporated into our simple fits. The dimuon with missing energy final 



state shows an excess of data compared to the standard model Monte Carlo prediction. A study 
into the track curvature of data and MC muons, and of the associated resolution, has shown 
that an additional smearing should be applied in the Monte Carlo to appropriately simulate 
very high pT muons. The prime signature of these muons is an excess of because of the lack 
of compensation for the mismeasured, unbalanced track. 

3 Targeted Searches using bump hunter (CDF) and sleuth 

The BUMP HUNTER algorithm scans the spectrum of most mass variables with a sliding window. 

The window varies across each final state and within final states themselves to account for 
changing detector resolution. For each combination of objects, the mass is calculated in each 
sliding window, and the probability to see a data excess as large as the one observed is noted. 
The window then slides by two times the mass resolution. In order for a bump to be defined, 
the region must contain at least five data events, and the side bands (regions on cither side of 
the region of interest and of half the width) must be less discrepant than the central region of 
interest. The smallest probability for each final state is determined and the overall significance 
of the most significant of these excesses is evaluated by performing pseudoexperiments. 

One bump is found to cross the discovery threshold. This is the mass of four jets in a final 
state with four jets oi^pT < 400. This discrepancy is attributed to the same three jet effect 
as was observed previously in the viSTA shape discrepancies. The same effect is observed in 
other final states of different jet multiplicities with a significance that falls below the threshold. 

In SLEUTH, ^PT in each channel is searched for a cut which maximizes the significance 
of the data excess over the SM backgrounds producing J2pt > cut. This significance is then 
corrected for the number of trials in both the number of possible positions of the cut in the 
histogram, and at a higher level, the number of final states Sleuth examines. The final corrected 
probability corresponds to the probability that an individual final state would produce one or 
more probabilities as small as observed. A significant output from Sleuth is defined as a state 
with corrected probability < 0.001 (equivalent to about 3 Gaussian standard deviations). 

CDF finds no discrepant state that crosses the discovery threshold. The most discrepant 
states at CDF contain same sign muon-elcctron pairs. The most significant of these causes a V 
of 0.085 which is still well above the threshold of 0.001. Sleuth by definition does not focus on 
the correlation among the most discrepant objects and instead focuses on just the single most 
discrepant state. 

The VISTA final states that show broad numerical excesses at D0 are found again with the 
SLEUTH algorithm, as would be expected. One additional distribution crosses the discovery 
threshold of "P < 0.001, where V is the probability after all trial factors: the final state that 
crosses the discovery threshold is + + ^t- Currently the evidence suggests that the 
muon tracking resolution is responsible for this discrepancy. A large fraction of the events in 
the tail of the Sleuth distribution have a muon with a very large pT and large missing energy. 
With the present modeling of muon resolution, straight track events are underrepresented in 
the standard model background estimation. This state has 46 data events in the tail of this 
distribution compared to only 17 predicted by the Monte Carlo. 

4 Conclusion 

In conclusion, CDF and D0 have performed broad searches for new physics using 2.0 fh~^ at 
CDF and 1.0 fh~^ at D0 . Hundreds of exclusive data final states and thousands of kinematic 
distributions were compared to the complete standard model background predictions at each 
experiment using the viSTA algorithm. No final states were discrepant at CDF and only four 
out of 180 exclusive final states showed a statistically significant discrepancy at D0 . Given 
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Figure 1: CDF and D0 distributions of final state counts and Kolmogorov-Smirnov probabilities converted to 



units of a before trials factor. The CDF final state and shape distributions are shown in Figs. 1(a) and 1(b) while 



the corresponding D0 distributions are shown in Figs. 1(c) and 1(d) 



the known modeling difficulties in all four final states, we refrain from attributing the observed 
discrepancies to new physics. A quasi-model independent search for new physics was also per- 
formed using the algorithm Sleuth by looking at the regions of excess on the high-^ tails of 
exclusive final states. Only /i^ -|- -|- Tpr surpasses the discovery threshold beyond the obvious 
excesses noticed in vista, and this seems to be related to difficulties in modeling the muon 
resolution. At CDF, there were no Sleuth discrepancies, but the most discrepant final states 
contained a same-sign muon and electron. Since no final state met the criteria defined before- 
hand for discovery, no new physics claim is made. Although there were no strong hints of new 
physics in the data, a factor of five more data has already been recorded at each experiment. As 
this data is incorporated into the analyses and improvements are implemented in the correction 
models, both experiments should be much more sensitive to possible new physics. 
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